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O \ Using the isospectral partners of the Bessel functions derived by Reyes et al. pQ, we 

find, on one hand, that these functions show non-typical supersymmetric (SUSY) 
behavior and, on the other, that the isospectral partner of the classical wave equation 
is equivalent to that of a damped system whose oscillations do not vanish in time, 
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but show a non-harmonic shape. 
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O r Usando las compaheras isoespectrales de las funciones de Bessel obtenidas por Reyes 

et al. p], encontramos, por un lado, que estas funciones muestran un compor- 
tamiento ati'pico de SUSY, mientras que, por otro lado, la compahera isoespectral 
de la ecuacion de onda clasica es equivalente a la de un sistema amortiguado cuyas 
oscilaciones no desvanecen con el tiempo, sino que obtienen una forma que no es 
armonica. 



Descriptores: Supersimetria; ecuacion de Bessel. 
^ ! PACS: 02.30.Hq; 03.65.-w; 11.30.Pb 



g ; 1 Introduction 

In quantum mechanics (QM), the number of problems which can be exactly solved is 
very limited, and one can only hope to get the approximated solution using a variety of 
methods, or turn to a modified problem which can be described in terms of the known 
exact problems. In fact, one of the main virtues of supersymmetric (SUSY) quantum 
mechanics is that there one can find an infinite number of one-parameter problems which 
posses the same spectra as the known exact ones. If SUSY cannot be applied, one can still 
find isospectral solutions by using the classical factorization method, or the interwinning 
method, to look for new equations with the old spectra. 

In classical mechanics, the only problem which can be directly compared to quantum 
mechanics is that described by the classical wave equation 

V 2 -\)^(x,y;t) = 0, (1) 
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but it has been stated before that there is no SUSY partner of this equation [2J. In fact, 
a factorization of the Bessel equation a la Infeld and Hull cannot be find [3]. However, 
the Bessel equation, 

d 2 JJr) IdJJr) / n T 




Jn(r) = 0, (2) 
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with n > 0, which arises from the wave equation after separation of variables, still possess 
a factorization in terms of raising and lowering operators defined by the equations [I] 

AUn(r) = - fj J n (r) = -J n+1 (r) , (3) 
A- +1 J n+ i(r) ={j r + J n+1 (r) = J n (r) , (4) 

respectively. 

Following the work of Mielnik [5j, who finds a family of potentials which posses the 
same spectrum as that of the one dimensional harmonic oscillator, and the work of Pina 
[6] on the factorization of some special functions found in mathematical physics, Reyes 
et al. [1] are able to derive second order differential equations which are 'isospectral' 
to the equations described in the Sturm-Liouville theory. Note that since the 'spectral' 
parameter n appears in and A~, contrary to SUSYQM [5], it shows up in the partner 
equation in a more complicated fashion than the original equation. 

As one can see, obtaining the isospectral partner of the Bessel equation is the first 
step toward an isospectral classical wave equation, this being the main purpose of this 
letter. But, before proceeding in this direction, we first show that the partner functions of 
the Bessel functions show a unique and unusual SUSY behavior. Then, we show that the 
isospectral classical wave equation that comes from the isospectral partner of the Bessel 
equation resembles the problem of damped waves, which nevertheless do not vanish in 
time, but show non-harmonic shapes due to this damping term. 



2 Isospectral Bessel Equation 

In ref.p], the isospectral partner of the Bessel equation was found as 

d 2 J n+1 l dJ n+1 ( (n + l) 2 \ ~ An {(2n + l) 7 r 2rt + 1] j 

dr 2 + r dr + { r 2 J n+1 ~ r 2 ( 7 r 2 « + l) 2 n+1 ' 1 J 

where, as is usual in SUSYQM, n > and the lowest lying eigenvalue is lost. Here, 7 is 
the parameter of the partner functions, with < 7 < 00. 

The term on the right of eq.([S]) corresponds to the extra potential function term in a 
typical SUSYQM problem, with the difference that in this case, the integral defining this 
term is exactly solvable. This also happens in the partner Bessel functions, which were 
found to be 

J n+1 (r; 7 ) = -J n+ i{r) + y (ir 2n + l)" 1 J n (r) . (6) 

Notice that the fact that these functions are explicitely found allows one to take a 
closer look to their properties, and by doing so we are here able to show that there exist 
a non-typical SUSY behavior for the isospectral Bessel functions, in the following sense. 
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In ref. [T] it was thought that the partner functions where not regular at r = 0. This is 
not the case, though. One can see that J\ = —J\, and that for n > we can use one of 
the recursion relations among Bessel functions, to say 

2n 

—J n (r) = J n+1 (r) + J n _i(r) , (7) 
r 

to find that the partner Bessel functions are regular at r = 0, since then 
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Now, this is a very unique feature of the Bessel partner functions. In SUSYQM the 
partner function of order n + 1, i/j n+ i(x; 7) , relate to the original function of preceding 
order ip n {x) through a differential operator. Here, the partner function of order n + 1 
relate to two Bessel functions, of orders n — 1 and n + 1 (see fig.([T]).) Moreover, one can 
see that the limits of the parameter 7 give 

Jn+i(r;7 = 0) = J„-i(r), (9) 

while 

J„+i(r;7 = oo) = -J n+ i(r) . (10) 

Therefore, the partner Bessel function J n+ i(r; 7), transforms from J n _i(r) to — J„+i( r ) 
as 7 goes from to 00. This is its most unusual characteristic, to be related to a Bessel 
function of two orders less than its own. 




Figure 1: Typical SUSY isospectral function generating scheme (left), and the way Bessel 
isospectral partners are generated (right), in terms of the n — 1 and n + 1 orders of the 
original functions. 
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3 Stationary oscillations in a damped wave equation 



In order to connect the partner Bessel funcions to the wave equation, we begin writing 
eq.(E]) for kr instead of r, 



2 d J dJi 



dr 2 



n+l 



dr 



where 



SWfi(u;7) 



+ (k 2 r 2 - (n + l) 2 ) J n+1 = g n +i{kr; 7) J n+1 
4n[(2n + l) 7 w 2ri + 1] 



(12) 
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Now, we multiply by a function of the polar angle 6, H{6), and assume that this 
function satisfies the equation 



d 2 H 



d6 

in order to write the wave equation 
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, +(n + l) 2 H = 



V 2 - 



ip n+1 (r,9;t;-f) 



(13) 



(14) 



where ^ n+1 (r, 6; t; 7) = J n+1 (kr;'y)H{e)e iut . 

This is the wave equation derived from the isospectral Bessel functions of ref . pQ , pos- 
sessing a damping term g n+ i(kr; 7)/r 2 which, though singular at r = 0, allows stationary 
solutions whose radial part is the partner Bessel function (jSJ). It is the kind of equation 
that may help on our comprehension of neuronal activity, where such equations appear 





Figure 2: Partner Bessel functions: J 2 (left), evolving from J , 7 = 0, to — J 2 , 7 = 00, 
and J3, evolving from J\ to — J3 (right). Note how the damping term affects the typical 
harmonic shapes of the Bessel functions, except when 7 = and 7 = 00. 

One may ask then, how does this damping term reflects in the stationary waves? The 
answer comes from our discussion above, about the way J n+ i is determined by J n -\ and 
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J n +i- Notice that for 7 = the damping term reduces to An/r 2 and is absorbed into the 
ordinary Bessel equation (J5J), and that for 7 = 00 the damping term reduces to zero. For 
other values of 7, the damping term modifies the Bessel functions, making them lose their 
typical harmonic shapes, as can be seen in fig. (2), where drastic changes in the harmonic 
shapes are seen as 7 increases towards infinity, specially for r close to zero. 

4 Conclusion 

In this letter we have shown that it is possible to find an isospectral partner of the classical 
wave equation. By using the isospectral partners of the Bessel functions from ref.p], we 
have shown that the arising wave equation contains a damping term, whose action does 
not destroy the stationary waves, but makes them acquire non-harmonic shapes. Since 
the Bessel equation can not be supersymmetrized [2], this may be the only way one can 
find an isospectral partner of the wave equation. 
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